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A comprehensive kinetic model of slow pyrolysis of biomass during a thermogravimetric analysis (TGA) 
has been developed, including the simulation of variable heating rates, composition estimation and struc¬ 
tural analysis of biomass. Biomass was assumed as a matrix of three solid global components (hemicel¬ 
lulose, cellulose and lignin) in which water and oil can be also present. 

Kinetics were based on an auto-catalytic model because it can simulate the degradation in cellulosic 

materials, as the cleavage of the biopolymers produce oligomers that accelerate the further depolymer¬ 
isation. The reaction pathway followed the Waterloo’s mechanism, which stablishes that all solid com¬ 
pounds decompose into volatiles and charcoal. This mechanism was completed by the vapourization 
of water and oil, and assuming that the formed charcoal can break into volatiles by a slow reaction. 
The set was solved by the 8th Runge-Kutta’s method and validated by the Simplex Nelder-Mead and 
Broyden-Fletcher-Goldfarb-Shanno’s methods. The development of this model has a high interest 
because it can help to understand how the conversion from biomass to biochemicals takes place. 

To assess what parameters can affect the thermal degradation of biomass pure polymeric samples of 
hemicellulose, cellulose and lignin and complex samples (seeds and woody biomass) were studied and 
fitted. Two types of operations were considered too. An isothermal degradation from 150 °C up to 
350 °C with increments of 50 °C, and a non-isothermal decomposition with heating rate of 5 °C/min, 
10 °C/min and 20 °C/min up to temperatures around 800 °C, depending on the studied sample. Average 
absolute deviations lower than 7% were obtained. It was deduced that there are some interactions 
between the three main biomass compounds. These interactions were observed by the variations in 
the kinetic parameters between complex and pure samples, also they were perceived between the iso¬ 
thermal and no-isothermal way. On the other hand, an effect of the biomass structure has been reported 
by the differences between the kinetics of the seeds and of the woody samples. It is remarkable that the 
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developed model could reproduce the cellulose decomposition with a variable heating rate using a 
unique set of kinetic parameters. This was possible by a no-Arrhenius’ dependence with temperature. 
In the same way, it was used to predict the initial composition of the studied biomass with deviations 
lower than 7% for lignin and cellulose. 

© 2015 Elsevier Ltd. All rights reserved. 


Nomenclature 



Acronyms 

h 

partial mass transfer coefficient between the liquid and 

C 

cellulose 


the gas, kgj m 3 /min m 2 kmoly 

HC 

hemicellulose 

k 0i 

preexponential factor for the reaction “i”, min -1 

L 

lignin 

ki 

kinetic constant for the reaction “i”, min -1 

O 

oil 

M e xp 

experimental mass fraction of unreacted biomass, g sam _ 

TGA 

thermogravimetric analysis 


ple/Ssample initial 

W 

water 

M 

mass fraction of unreacted biomass, g S ampie/gsam P ie initial 

XC 

charcoal produced form cellulose 

m Jm ax 

maximum value for mass fraction of the compound “j” 

XHC 

charcoal produced form hemicellulose 


in the biomass, g/g 

XL 

charcoal produced form lignin 

m j m in 

minimum value for mass fraction of the compound “j” 




in the biomass, g/g 

Subindex and sunerindex 

m i 

mass fraction of the compound “j” in the biomass, g/g 

EXP 

experimental data of the TGA 

N 

number of compounds in the biomass, adim 

in 

inert compounds 

n i 

Order of reaction of the reaction “i”, adim 

TOTAL 

total simulated TGA 

nli 

mass transfer order, adim 



N r 

number of reactions, adim q 

Greek letters and symbols 

?! 

vapour pressure of the compound “j”, atm 

OCi 

initialization factor, adim 

n 

reaction velocity number “i”, g/min g 

Pi 

acceleration factor, adim 

0 

reaction velocity of decomposition for the component 

rl 

Ai - Fj 

Antoine’s equation coefficients of the compound “j”, 


“j” in the biomass, g/min g 

j j 

adim 

s 

exchange surface between the liquid and the gas, m 2 

c 

correction factor for the kinetic in the decomposition at 

t 

operating time, min 


different heating rates of the cellulose, adim 

T 

operating temperature, K 

C 

concentration of “j” in the gas phase, kmol/m 3 

*iexp 

experimental biomass fraction, g S ampie/g S am P ie initial 

j 

c; 

equilibrium concentration of “j” in the interphase 

X isiM 

simulated biomass fraction, g sam pie/gsampie initial 


between the liquid and the gas phase, kmol/m 3 



Ea t 

~R 

activation energy of the reaction “i”, K 




1. Introduction 

The use of fossil fuels as the main raw material for industry is not 
sustainable, and certainly it will not be the forever-solution. So a 
new source of basic compounds (i.e. carbon, hydrogen and oxygen) 
and energy should be considered. This new source could be biomass 
[1 ], which can be transformed into bioenergy, biochemical and bio¬ 
fuels in biorefineries [2,3]. However, the design of these biorefiner¬ 
ies requires knowledge about the conversion from raw material to 
fuels and fast, cheap and accurate biomass-analysing methods. For 
the latter, several wet methods of chemical analysis have been used 
[4]. These methods are based on the fractionation of biomass sam¬ 
ples and a later isolation of purified fractions, which could be quan¬ 
tified using conventional analytical instruments. Although these 
techniques have high accuracy and robustness, they are not suitable 
for an industrial scale because they are expensive and require a lot 
of time. Another option would be spectroscopic analysis, such as, 
the Near Infrared Reflectance (NIR) spectroscopy, which reduces 
time requirements and cost and it is a method with a high repro¬ 
ducibility. Nevertheless, these analysis need data with a very high 
quality and an initial blank spectrum, which is an important limita¬ 
tion. So, the measurement of the initial biomass composition is an 
issue that have not an optimal solution yet. Thermogravimetric 
analysis (TGA) of biomass could be the answer for this problem 
under certain conditions. In addition, it can provide information 
about how the thermal decomposition takes place. 


Thermogravimetric analysis is a slow pyrolysis process which 
consist of recording the mass variation of a sample which is treated 
with a temperature profile. This profile is provided by a gas phase 
which can be an inert or an oxidant compound [5]. This type of 
analysis have been studied thoroughly [5-7] and there are a lot 
of works about their modelling in the literature. The most 
extended model considers a first order kinetic for each compound 
present in biomass assuming that biomass is formed by three main 
compounds (cellulose, hemicellulose and lignin). These compo¬ 
nents decompose to charcoal and volatiles by independent reac¬ 
tions. Volker and Rieckmann [8] used a first order kinetic to 
adjust the decomposition of pure cellulose and the deviation 
between the experimental data and the simulation was relatively 
high. In contrast, Capart et al. [9] studied the pure cellulose ther¬ 
mal breaking but considering an autocatalytic model which sup¬ 
posed a good fitting with an overall deviation around 1%. On the 
other hand, Mangut et al. [10] proposed a kinetic model of n th - 
order for the degradation of residues from tomato processing 
industry which could reproduce the biomass behaviour. But Zaba- 
niotoua et al. [11], Slopiecka et al. [12] and Kastanaki et al. [13] 
studied the TGA kinetics of several lignocellulosic biomass sam¬ 
ples, poplar wood and lignite-biomass blends respectively with a 
first order model and they obtained good fits too. Slopiecka et al. 
and Zabaniotoua et al. fitted their TGA as a single compound, which 
is useful to reproduce the decomposition. However, it is not 
capable of reproducing the individual behaviour of the biomass 
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Table 1 

Initial composition of the samples. 



Water 

(g/g) 

Oil 

(g/g) 

Hemicellulose 

(g/g) 

Cellulose 

(g/g) 

Lignin 

(g/g) 

CharHemicellulose 

(g/g) 

Char Ce iiuiose 

(g/g) 

Char Lignin 

(g/g) 

Inert 

(g/g) 

Hemicellulose 

0.030 

0.000 

0.970 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Cellulose 

0.009 

0.000 

0.000 

0.991 

0.000 

0.000 

0.000 

0.000 

0.000 

Lignin (Turku) 

0.011 

0.000 

0.000 

0.000 

0.989 

0.000 

0.000 

0.000 

0.000 

Lignin (Alkaline) 

0.038 

0.000 

0.000 

0.000 

0.962 

0.000 

0.000 

0.000 

0.000 

No extracted seeds 

0.029 

0.166 

0.146 

0.214 

0.419 

0.000 

0.000 

0.000 

0.026 

Extracted seeds 

0.015 

0.000 

0.180 

0.264 

0.515 

0.000 

0.000 

0.000 

0.027 

Hydrolysed seeds 

250 °C 

0.019 

0.000 

0.054 

0.239 

0.463 

0.104 

0.030 

0.050 

0.041 

300 °C 

0.019 

0.000 

0.000 

0.259 

0.463 

0.110 

0.040 

0.069 

0.041 

340 °C 

0.019 

0.000 

0.000 

0.029 

0.416 

0.110 

0.248 

0.090 

0.088 

Seed skin 

0.014 

0.100 

0.360 

0.439 

0.000 

0.000 

0.000 

0.000 

0.087 

Lime 

0.014 

0.049 

0.215 

0.494 

0.209 

0.000 

0.000 

0.000 

0.019 

Plane tree 

0.014 

0.005 

0.220 

0.501 

0.260 

0.000 

0.000 

0.000 

0.000 

Holm oak 

0.017 

0.023 

0.196 

0.455 

0.308 

0.000 

0.000 

0.000 

0.001 

Cedar 

0.021 

0.035 

0.230 

0.331 

0.382 

0.000 

0.000 

0.000 

0.000 

Catalpa 

0.014 

0.009 

0.230 

0.468 

0.261 

0.000 

0.000 

0.000 

0.018 

Acer Saccharum 

0.015 

0.015 

0.213 

0.426 

0.329 

0.000 

0.000 

0.000 

0.000 

Almond tree 

0.018 

0.015 

0.200 

0.400 

0.366 

0.000 

0.000 

0.000 

0.002 


components and ensuring that the obtained parameters have phys¬ 
ical meaning. On the other hand, Kastanaki et al. and Mangut et al. 
adjusted their TGA with individual kinetics for each biomass com¬ 
pound. Therefore, the behaviour of each of them could be simu¬ 
lated and the physical sense of the parameters checked. 
Nevertheless, they did not studied the causes of the variations in 
the kinetics of the biomass components assuming that there are 
no interactions between them. Taking into account this big range 
of possible models it is difficult to select one because any of them 
could be a good way to simulate the thermal degradation of bio¬ 
mass. Finally, the autocatalytic model is the option selected in this 
work due to the fact that it can reproduce the steep changes in cel- 
lulosic material better than a first or nth-order model. 

The aim of this paper is to study the main parameters that affect 
to the thermal decomposition of lignocellulosic biomass in a 
thermogravimetric analysis developing a model which could 
reproduce the decomposition of any sample of biomass and its 
components in an isothermal or non-isothermal process and at 
any heating rate. The last one is important because a discrepancy 
in the amount of formed charcoal has been reported by other 
authors [8,9,12,14,15] when different heating rates are used. Fur¬ 
thermore, the causes of the changes in the kinetics of the individ¬ 
ual biomass components with the type of samples or process was 
assessed by the comparison of the kinetics parameters, which is 
not generally done in previous studies. In addition, this model 
should be able to estimate the initial composition of the degraded 
sample from the thermogravimetric analysis and from the kinetic 
parameters fitted previously. This capacity is important because 
it is a new use for the TGA modelling and, if it is developed cor¬ 
rectly, it would become an economic option to obtaining the initial 
composition of the biomass. 




Fig. 1 . Schema of the biomass structure. 


2. Materials and methods 

2.1. Materials 

Grape seeds from Vitis vinifera L (Tempranillo) from Matarro- 
mera S.A. winery (Valbuena de Duero, Spain) campaign 2011 and 
several woody wastes were used as raw material. Material for this 
study was ground to a particle size of 0.5-1.0 mm. 

The reagents used for HPLC analysis were: cellobiose (+98%), 
glucose (+99%), fructose (+99%), glyceraldehyde (95%), pyruvalde- 
hyde (40%), arabinose (+99%), 5-hydroxymethylfurfural (99%), 
lactic acid (85%), formic acid (98%), acrylic acid (99%), mannose 
(+99%), xylose (+99%) and galactose (+99%) purchased from Sigma 
and used without further modification. For the structural 
carbohydrates and lignin determination sulphuric acid (98%) and 
calcium carbonate (^99.0%) were used as reagents supplied by 
Sigma too. 


2.2. Biomass characterization 
2.2A. Sugar content 

The sugar content measurement of a biomass sample requires 
its hydrothermal fractionation followed by a hydrolysis of the 
product which will be fed to a HPLC later. The hydrolysis is needed 
because the fractionation generates a range of polymeric fractions 
and compounds which could not be directly identified in a HPLC. 
So they have to be cleavaged by a further hydrolysis into their basic 
units or monomers, e.g. glucose, fructose, xylose and arabinose. 

The samples were hydrolysed adding 3.00 ± 0.01 ml of sulphu¬ 
ric acid (72%) to 15 ml of each aliquot. Each sample was incubated 
in a forced convection oven for 30 ± 5 min at 30 ± 3 °C. After this 
time, the samples were taken out from the oven, they were diluted 
by 84.00 ± 0.04 ml of deionized water and finally they were placed 
in the oven for 1 h at 121 °C. Afterwards, the solution was cooled 
down to room temperature and it was filtered under vacuum. 
Before injecting in the HPLC the samples were neutralized to 
pH = 6-7 using calcium carbonate. 

The HPLC column used for the separation of the compounds was 
a SUGAR SH-1011 Shodex at 50 °C at a flow of 0.8 ml/min using a 
solution of 0.01 N of sulphuric acid and Milli-Q water as mobile 
phase. A Waters IR detector 2414 and Waters dual X absorbance 
detector 2487 (210 nm and 254 nm) were used to identify the sug¬ 
ars and their derivatives. 
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Fig. 2. Reaction pathway in a thermal decomposition. 


2.2.2. Solid analysis. Klason lignin determination and sugars attached 
to the solid 

The raw material and the solid residue generated by the hydro¬ 
lysis were analysed for lignin content using the Klason assay 
according to the TAPPI standard method T-222 om-98 [16]. To 
do so, 300 mg of sample was put into laboratory glass bottles, 
3 mL of sulphuric acid (72%) was added and it was incubated dur¬ 
ing 30 min at 30 °C and it was shaken vigorously every 5-10 min. 
Then, the mixture was diluted with 84 mL of deionized water and 
it was placed in an oven for 1 h at 121 °C. At that moment, the 
sample was taken out from the oven, cooled down to room tem¬ 
perature and the mixture was filtered under vacuum. The 
obtained solid after filtration was dried at 105 °C for 24 h, it 
was cooled down in a desiccator and then it was weighted. This 
solid was introduced in the calcination oven at 550 °C for 24 h 
to determine the ash content. Considering the weight differences, 
the Klason lignin content was calculated. The hydrolysis liquid 
was neutralized with calcium carbonate to pH = 6-7, then it was 
filtered and analysed by HPLC as explained in Section 2.3.1.2 
Sugars. 

The initial composition calculated by the methods described in 
2.2.1 and 2.2.2 is collected in Table 1. 

2.3. Experimental set-up and procedure 

TGA were carried out in a TGA/SDTA RSI analyzer of Mettler 
Toledo. Samples of approximately 10 mg were heated from 50 °C 
to the required temperature at a rate of 20 °C/min under N 2 atmo¬ 
sphere (60 N mL/min flow) to determine the carbonization. The 
final temperature changed with the type of analysis. If the study 
was at isothermal conditions it had a value between 150°C and 
350 °C. However, when it was no isothermal the biomass were 
heated up to temperatures around 800 °C. 

2.3.1. Procedure for the analysis of the effect of the composition 

Thermogravimetric analysis at variable temperature with a 
heating rate of 20 °C/min of woody samples with different lignin 
content were fitted. So the difference between the kinetics param¬ 
eters were used to discover how the composition affects to the 
thermal degradation. TGA of hemicellulose, cellulose and grape 
skins, which does not have lignin, at the same heating rate were 
performed to study this factor too. 


Table 2 

Antoine’s equation coefficients of water (W) and oil (O). 



W 

O 

A 

6.21 E+01 

1.22E+01 

Bj 

—7.26E+03 

5.88E+03 

Q 

0.00E+00 

-2.93E+02 

Dj 

7.30E+00 

0.00E+00 

Bj 

4.17E-06 

0.00E+00 

Bj 

1.00E-02 

0.00E+00 


2.3.2. Procedure for the analysis of the effect of the structure 

The effect of the biomass structure was studied comparing the 
adjusted kinetic parameters obtained from the TGA (with a heating 
rate of 20 °C/min) of two types of pure lignin: an alkaline lignin and 
a sample from Turku, Finland. The last one was extracted using a 
hydrotropic substance, the p-toluene sulfonate. In addition, the 
deviation of the kinetics parameters between the TGA (heating rate 
of 20 °C/min) of a sample of grape seeds and grape seeds extracted 
with a mixture of ethanol/water (70/30) for 1 h was considered. 
The kinetics variation between these grape seeds and hydrolysed 
grape seeds for 1 h at three different temperatures (250 °C, 
300 °C and 340 °C) and at a heating rate of 20 °C/min were ana¬ 
lysed too. 

2.3.3. Procedure for the analysis of the effect of the heating rate 
The role of the heating rate was considered by fitting TGA of 

pure cellulose at three different heating rate: 5 °C, 10 °C, 20 °C 
and comparing the values of their kinetics parameters. 

2.3.4. Procedure for the analysis of the effect of the isothermal 
conditions 

This factor was studied by the adjustment of the TGA of Acer 
Saccharum, a type of maple, at 5 temperatures (150 °C, 200 °C, 
250 °C, 300 °C and 350 °C) with a heating rate of 20 °C/min and 
considering the modifications in the kinetics. 

3. Mathematical model 

3.1. Biomass composition 

The solid organic compounds present in biomass are divided in 
three main biopolymer fractions: hemicellulose, cellulose and 
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lignin. Hemicellulose and cellulose are constituted by monomeric 
sugars. The difference between them is that cellulose is a linear 
polymer of anhydroglucopyranose (hexose) units linked by ether 
bonds, while hemicellulose is a branched and amorphous polymer 
formed by both pentoses and hexoses. In contrast, lignin is a com¬ 
plex, cross-linked, three-dimensional aromatic polymer formed 
with phenylpropane units [17]. The lower the amount of lignin 
the more flexible is the biomass (e.g. herbs). In addition, some 
inerts (inorganic mineral compounds) and two liquid phases can 
be present inside the biomass: water and an organic phase which 
is identified in this work as oil. 

Simplifying the raw biomass structure we can consider that the 
cellulose microfibers are connected by hemicellulose in 3D struc¬ 
ture of lignin that encloses and protects them (Fig. 1). The structure 
of the biopolymer fractions and other kinds biomass lignin-lean 
can be a bit different. 


3.2. Reaction pathway 
3.2A. Degradation process 

The thermal degradation of biomass in an inert atmosphere 
(slow pyrolysis conditions) starts with the vapourization of liquid 
phases. Water evaporates near 100°C and oil between 100 °C 
and 300 °C. In this research, we intentionally did not dry the bio¬ 
mass until full dryness to mimic as much as possible some kind 
of humid conditions and therefore the water evaporation. 

At 200 °C lignin begins its decomposition, breaking its weaker 
parts and enhancing the reaction of hemicellulose and cellulose. 
Between 250 °C and 275 °C hemicellulose reacts and around 
300 °C it disappears completely, which promotes the cellulose 
breaking. The last one commences its degradation between 
300 °C and 350 °C and, from this point, only lignin, inert substances 
and the product from the decomposition (charcoal) remains in the 
biomass. Lignin depletes at 500 °C and charcoal continues in the 
sample with a very low degradation rate. Charcoal only fades com¬ 
pletely if the atmosphere is changed to an oxidant compound. 

3.2.2. Reaction mechanism 

The most extended idea in the literature [5,6,9-13] is a reaction 
pathway based on the Waterloo-mechanism. This theory estab¬ 
lishes that, under slow pyrolysis conditions, biomass decomposes 
to charcoal and gases via dehydration reactions [14,15]. Following 
this theory, it was assumed in this work that any organic solid 
compound in the biomass could be converted into charcoal and 
volatiles by independent parallel reactions. Furthermore, an indi¬ 
vidual charcoal for each solid organic component was used in this 
paper (Fig. 2), with the intention of describing the slow pyrolysis 
process in detail. This leads to a substantial improvement in the fit¬ 
ting of the experimental data better further prediction. The path¬ 
way was completed adding the decomposition of each charcoal 
to volatiles and with the vapourization of liquid phases if present. 


Table 3 

Averaged mass transfer parameters of water (W) and oil (O). 



h (kg m 3 /min m 2 kmol) 

nli 

Pure 

W 

1000 

2.00 

O 

- 

- 

Grapes 

W 

1000 

2.00 

O 

123 

1.00 

Wood 

W 

3000 

2.00 

O 

123 

1.00 


Table 4 

Groups of woody samples taking into account its lignin content. 


Lignin content (wt%) 

Samples 

A.A.D. a (%) 

20 

Lime 

1.97 

26 

Plane tree 

1.46 


Catalpa 

1.87 

30 

Holm oak 

1.82 

>30 

Cedar 

1.52 


Acer Saccharum 

2.14 


Almond 

4.31 

a Average absolute deviation between experimental and simulated data. 


Table 5 

Initial composition variation ranges for the composition estimation of no-extracted 
grape seeds. 


m m in a (§/§) 

m max b (g/g) 

W c 

0.00 

0.08 

O d 

0.00 

0.20 

HC e 

0.10 

0.25 

C f 

0.15 

0.60 

L s 

0.15 

0.45 


a The lowest mass fraction in the optimization. 
b The highest mass fraction in the optimization. 
c Water content. 
d Oil content. 
e Hemicellulose content. 
f Cellulose content. 
g Lignin content. 


3.3. The model 
3.3A. Assumptions 

Aimed at simplifying the modelling problem it was assumed 
that: 


a. All the reactions are irreversible and independent. So, the 
degradation kinetics of each component only depend on 
their composition and temperature [10,12,13]. 

b. There are no energy transport limitations within the bio¬ 
mass particles (as only 10 mg of micronized particles were 
used for the TGA analysis). Consequently, all the parts of 
the biomass are at the same temperature. 

c. Diffusional mass transport resistances for liquid phases are 
negligible, as the particles were micronized. 


3.3.2. Mass balances 

The model of the decomposition process considered a non-sta- 
tionary mass balance for each component in the biomass sample: 


dm, 

~dt~ Vi ~ 



( 1 ) 


And the total variation of mass was calculated by the addition of 
all of them: 

dM_J^dmj 

~dt~j^~dt ( j 

3.3.3. Kinetics 

Mass variation was caused by two different phenomena: reac¬ 
tion kinetics for the solid material and mass transfer of the liquid 
phases. 


• Liquid mass transfer: free oil (or extractable compounds) and 
free water. 
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Table 6 

Comparison between the estimated and experimental composition of no-extracted grape seeds. 



W a 

O b 

HC C 

C d 

L e 

in f 

m (wt%) g 

Experimental 

0.0292 

0.1655 

0.1461 

0.2142 

0.4187 

0.0263 

Estimated 

0.0469 

0.1149 

0.1887 

0.2215 

0.4121 

0.0159 

Deviation (%) h 

60.5% 

-30.6% 

29.2% 

3.39% 

-1.57% 

-39.5% 


a Water. 
b Oil. 

c Hemicellulose. 
d Cellulose. 
e Lignin. 
f Inert. 

s Biomass composition in weight percentage. 
h Deviation between the estimated and real composition. 


Mass transfer of liquid phases was described by the partial mass 
transfer coefficient in the gas phase, the mass transfer area and the 
difference between de equilibrium concentration in the liquid 
phase and the global concentration in the gas phase (as driving 
force) (3). 

r t = h -S- (C; - Cj) (3) 

As the operating pressure is the atmospheric, the equilibrium 
concentration (4) was obtained by the ideal gas equation and the 
vapour pressure calculated by the Antoine equation (5) of each 
compound. The Antoine’s equation coefficients of each liquid phase 
are compiled in Table 2. 


In (PJ) = A,- + JL + Dj ■ In(T) + Ej ■ T F > (5) 

In addition, the transfer area was considered as a function of the 
mass in the solid, so the final expression for the mass transfer was: 

r i = h-{q)-mf (6) 

• Solid kinetics: for solid organic compounds 
As it was mentioned in the introductory section, there are two 
options for temperature dependent kinetics, i.e. a first order reac¬ 
tion (7) and an autocatalytic reaction (8). The first is the most 
extended option in the bibliography [9-13] and the second is pro¬ 
posed because its response is very similar to the behaviour of the 
biomass observed in the literature [9] and in previous studies. Both 
kinetic equations considered an Arrhenius’ dependence with 
temperature. 

Ea t 

h = ki • irii = ko t • e~*T ■ m, (7) 


sition of cellulose was simulated by equation (10). This other kind 
of reaction rate expression was needed because the biomass shows 
a different behaviour when the heating rate varies along with time 
(or temperature), as the polymeric structure can both collapse or 
swell [8,9,14,15]. The parameter “c” is a correction factor to the 
cellulose decomposition when a variable heating rate is used. Nev¬ 
ertheless, we assumed that it could be affected by the biomass 
structure and heating process too. For this reason, three sets of 
“c” values are present in Table 7. 

r, = k 0i • e _ w +cT+ln(T) • mj 1 ' • (1 - a,- • mj)^ (9) 


3.4. Resolution 

The system of 8 ordinary differential equations (ODE) that 
results from the model was solved by the Runge-Kutta’s method 
with 8th order of convergence [18]. The validation of the model 
with the experimental data was done applying the Simplex Nel- 
der-Mead method for obtaining an initial estimation of the param¬ 
eters together with the Broyden-Fletcher-Goldfarb-Shanno’s 
method to improve this initial solution [18]. During the optimiza¬ 
tion, the optimization range were selected in order to achieve an 
optimum which would have physical meaning. 

The objective function used is the Absolute Average Deviation 
(AAD) defined as follows: 


aad = £ 


I^i'exp % i s 


^•100 


( 10 ) 


The developed program is available for free in the web page of 
the research group of high pressure processes of the University of 
Valladolid (http://hpp.uva.es/software/). 


4. Results and discussion 


n = k 0i ■ e~^r ■ m" 1 ■ (1 - cci ■ mjf (8) 

The coefficient a,- is the initialization factor which indicates the 
resistance of the biomass against the degradation. It is used to 
establish the initial value of the reaction velocity. In this work a, 
was fixed at 0.99, as it is the most recurrent value in the literature 
[9]. The coefficient is the acceleration factor and represents how 
fast the degradation is once it has started. The autocatalytic kinet¬ 
ics can predict the dramatic changes in the total mass along with 
temperature better than a first order kinetics. In view of that, auto¬ 
catalytic kinetics was the selected option for this work. 

In addition, as pure cellulose has been studied at different rates 
of heating a non-Arrhenius’ dependence with the temperature was 
added to equation (9) for this compound. Therefore, the decompo- 


4A. Pure samples 

This adjustment was done taking into account the theoretical 
development showed in Section 3. The initial composition of each 
compound is arrayed in Table 1. 

4AA. Hemicellulose 

The fitting of pure hemicellulose decomposition (with an aver¬ 
age absolute deviation of 2.2%) is shown in Fig. 3 and the kinetic 
parameters in Table 7. The mass transfer parameters were aver¬ 
aged for all the samples (complex samples too) and are shown in 

Table 3. 

It can be observed in Fig. 3 that there was a first reduction 
between 50 °C and 105 °C which corresponds to water vapouriza- 


Please cite this article in press as: Cabeza A et al. Autocatalytic kinetic model for thermogravimetric analysis and composition estimation of biomass and 
polymeric fractions. Fuel (2015), http://dx.doi.Org/10.1016/j.fuel.2015.01.048 
















Please cite this article in press as: Cabeza A et al. Autocatalytic kinetic model for thermogravimetric analysis and composition estimation of biomass and 
polymeric fractions. Fuel (2015), http://dx.doi.Org/10.1016/j.fuel.2015.01.048 


Table 7 

Kinetics parameters fitted for all the samples. 



koi 

ko 2 

ko 3 

ko 4 

ko 5 

ko 6 

ko 7 

ko 8 

kOg 

Ea^R Ea 2 /R 

Ea 3 /R Ea 4 /R Ea 3 /R Eag/R Ea 7 /R Ea 8 /R Eag/R Pi 

h 

(h 

Pa 

P 5 

@6 

p 7 

(h 

Ps 

c 

Hemicellulose 

19,894 

_ 

_ 

_ 

_ 

_ 

20,962 

_ 

0.024 

7,177 

_ 

- 

_ 

_ 

_ 

7,120 

_ 

101 

0.0000 

- 

- 

- 

- 

- 

0.7700 

- 

0.5200 

_ 

Cellulose 

- 

5,976 

- 

- 

- 

32,613 

- 

0.014 

- 

- 

12,190 

- 

- 

- 

8,658 

- 

274 

- 

- 

0.5300 

- 

- 

- 

0.0000 

- 

1.1300 

- 

0.0060 

Lignin 

- 

- 

48,433 

48,857 

0.013 

- 

- 

- 

- 

- 

- 

8,308 

8,192 

481 

- 

- 

- 

- 

- 

- 

0.0000 

0.0000 

0.0000 

- 

- 

- 

- 

- 

(Turku) 

Lignin 



33,297 

14,151 

0.012 







7,579 

6,482 

122 







0.0000 

0.0000 

0.0000 






(Alkaline) 





























No extracted 
seeds 

Extracted 

seeds 

Hydrolysed 

seeds 

Grape skin 

19,894 

5,976 

48,594 

47,982 

0.015 

32,613 

20,962 

0.014 

0.024 

7,177 

12,190 

9,337 

8,936 

100 

8,658 

7,120 

274 

101 

0.0000 

0.5300 

0.0000 

0.1100 

1.9200 

0.0000 

0.7700 

1.1300 

0.5200 

0.0060 

19,580 

6,023 

48,619 

48,046 

0.045 

32,629 

22,726 

0.014 

0.024 

6,097 

12,120 

8,677 

8,438 

100 

8,570 

6,009 

274 

101 

0.0050 

0.5300 

0.0000 

0.1000 

1.9100 

0.0000 

0.7600 

1.1300 

0.5200 

0.0060 

19,580 

5,648 

12,148 

48,415 

0.002 

32,480 

22,726 

0.014 

0.024 

6,097 

13,286 

9,625 

9,898 

124 

9,248 

6,009 

274 

101 

0.0050 

0.5300 

0.0000 

0.1000 

1.9100 

0.0000 

0.7600 

1.1300 

0.5200 

0.0060 

17,621 

4,221 

_ 

_ 

_ 

40,014 

21,272 

0.021 

0.067 

6,049 

12,269 

_ 

_ 

_ 

7,505 

6,584 

104 

103 

0.0050 

0.0714 

_ 

_ 

_ 

0.1021 

0.0152 

0.0079 

0.0134 

0.0040 

Trees-20% 

39,194 44,605 

48,812 

48,191 

0.018 

30,071 

7,120 

0.078 

0.024 

6,755 

11,127 

8,503 

8,104 

100 

7,798 

7,320 

274 

101 

0.6584 

0.7149 

0.0017 

0.1068 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0001 

Trees-26% 

39,194 43,605 

48,812 

48,191 

0.018 

30,071 

7,120 

0.078 

0.024 

6,755 

11,127 

8,203 

8,004 

100 

7,998 

7,320 

274 

101 

0.6584 

0.7149 

0.0017 

0.1068 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0001 

Tress-30% 

39,194 38,595 

47,812 

48,202 

0.018 

30,071 

7,120 

0.078 

0.024 

6,755 

11,105 

8,183 

7,602 

100 

8,098 

7,320 

274 

101 

0.6584 

0.7163 

0.0017 

0.1062 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0001 

Trees->30% 

39,194 27,580 

46,988 

46,970 

0.036 

29,653 

7,120 

0.082 

0.024 

6,755 

10,934 

8,870 

8,227 

100 

9,162 

7,320 

274 

101 

0.6584 

0.7163 

0.0017 

0.1062 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0001 

Isothermal 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

150-200°C 

17,028 

38,652 

8,120 

6,016 

0.000 

33,233 

43,823 

0.013 

0.000 

7,295 

13,210 

8,770 

8,060 

121 

8,544 

7,866 

390 

100 

0.6584 

0.7149 

0.0017 

0.1068 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0037 

250 °C 

16,973 

38,652 

8,120 

6,016 

0.000 

33,233 

43,825 

0.013 

0.004 

7,295 

13,210 

8,770 

8,060 

121 

8,544 

7,866 

390 

100 

0.6584 

0.7149 

0.0017 

0.1068 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0037 

300 °C 

16,525 

38,652 

8,120 

6,016 

0.001 

33,233 

43,825 

0.013 

0.918 

6,597 

13,210 

8,770 

8,060 

121 

8,544 

6,847 

390 

100 

0.6584 

0.7149 

0.0017 

0.1068 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0037 

350 °C 

16,525 

34,637 

10,144 

6,155 

0.001 

40,573 

43,825 

0.013 

0.918 

6,229 

13,210 

8,770 

8,060 

121 

7,544 

6,847 

390 

100 

0.6584 

0.7149 

0.0017 

0.1068 

1.9150 

0.0000 

0.7680 

1.1309 

0.5192 

0.0037 
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T,2C 

Fig. 5. Simulated differential thermography of the lime TGA. W: water. HC: 
hemicellulose. XHC: char of hemicellulose. TOTAL: simulated DTG. O: oil. C: 
cellulose. L: lignin. XC: char of cellulose. XL: char of lignin. 



Fig. 6. Simulated differential thermography and experimental differential ther¬ 
mography. EXP: experimental DTG. TOTAL: simulated DTG. 


example, Elyounssi et al. [14] and Slopiecka et al. [12] found that 
this peak is between 260 and 299 °C at low heating rates and Wil¬ 
liams and Besler [5] discovered that it is around 310 °C at 20 K/ 
min. This decoupling could be caused by the fact that the value 
obtained in this work is for extracted hemicellulose and the others 
from complex samples. However, the temperature of the hemicel¬ 
lulose maximum in our woody biomass fittings (304 °C) agrees 
with the value of these authors (Fig. 5). So, this deviation would 
show that the same compound in a different structure has a differ¬ 
ent behaviour against thermal degradation. Finally, the comparison 
between experimental and simulated DTG were done (Fig. 6). Both 
curves presents the same behaviour, with the main variations at 
the same temperatures. The differences between them could come 
from the numerical evaluation of the experimental derivative, 
which is obtained by a central difference approximation. 


4A.2. Cellulose 

The adjustment of the TGA of cellulose at heating rates of 20 °C/ 
min (average absolute deviation of 6.6%), 10 °C/min (average abso¬ 
lute deviation of 3.7%) and 5 °C/min (average absolute deviation of 
1.2%) was done in this point. The obtained kinetic parameters are 
presented in Table 7. The thermolysis process started at 300 °C, 
with the fractionation of cellulose, and continued until 400 °C. 
After this temperature, there was only charcoal in the sample. This 
experimental behaviour of cellulose agrees with the work of previ¬ 
ous authors [4,5,7,9,11,12] and it could be simulated using the 
same set of kinetics parameter for the three experiments. 


A decoupling was observed in mass variation originated by the 
different heating rates: the lower the heating rate is, the higher 
yield of charcoal is obtained. The origin of this fluctuation could 
be that a slow heating rate provides enough time to appear second¬ 
ary reactions which increase the charcoal production [9,14,15] due 
the collapse of the polymeric structure. For this reason, a specified 
model was proposed to cellulose decomposition (point 3.3). The 
idea was to represent this change in the decomposition process 
by a non-Arrhenius’ dependence with temperature. 

It is observed by the differential thermography of cellulose that 
its pure fraction shows a peak at 370 °C. Again, it decreased to 
346 °C in a woody sample (Fig. 5). This last value is similar to the 
data provided by Williams and Besler [5] and Carrier et al. [4]. Con¬ 
sequently, there is a modification in the thermal degradation of the 
cellulose when it is inside a complex sample due to the variation in 
structure. 

4A.3. Lignin 

Both studied samples showed a similar behaviour that the pre¬ 
viously described in the literature [4,5,7,11,12], the reaction began 
around 200 °C with the cleavage of the lignin and went on until 
500 °C when only charcoal remained in biomass. However, both 
lignin samples have some quantitative differences which did not 
allow to fit the experimental data with only one set of kinetic 
parameters. This discrepancy between them implies that it is very 
important how the biomass has been treated and how its structure 
is. The kinetic parameters from the adjustment of the lining from 
Turku (average absolute deviation of 1.7%) and from the fitting of 
the alkaline lignin (average absolute deviation of 1.1%) are listed 
in Table 7. 

The differential thermogravimetric analysis of the lignin 
showed a minimum at 322 °C for its isolated fraction. This value 
was far from the value of 387 °C, that appeared in a woody sample 
(Fig. 5) and from the value established by Williams and Besler [5] 
(390 °C) or Carrier et al. [4] (395 °C). This decoupling would be 
again caused by an effect of the structure of the sample. 

4.2. Complex samples 

In this section the TGA of complex samples was studied. The 
procedure was to try to simulate the experimental data from the 
TGA of complex samples with the kinetics parameters obtained 
for the pure fractions. Hence, only if this simulation had no corre¬ 
lation with the experimental performance, a new set of parameters 
would be searched. 



Fig. 7. Variation of the kinetic constants between the extracted seeds and the no 
extracted seeds. Kl: kinetic constant of hemicellulose degradation to volatiles. K2: 
kinetic constant of cellulose degradation to volatiles. K3: kinetic constant of lignin 
degradation to volatiles. K4: kinetic constant of lignin degradation to char. K5: 
kinetic constant of lignin char degradation to volatiles. I<6: kinetic constant of 
cellulose degradation to char. K7: kinetic constant of hemicellulose degradation to 
char. 
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542 4.2A. Grape seeds 

543 This section is focused in the TGA of wastes from the wine 

544 industry: seeds, extracted seeds, seed skins and hydrolysed seeds. 

545 4.2A A. No-extracted seeds and extracted seeds. The TGA of the no- 

546 extracted seeds was fitted (average absolute deviation of 1.6%) 



Fig. 8. Variation of the kinetic constantan between the hydrolysed seeds and the 
non-hydrolysed seeds. Kl: kinetic constant of hemicellulose degradation to 
volatiles. K2: kinetic constant of cellulose degradation to volatiles. K3: kinetic 
constant of lignin degradation to volatiles. K4: kinetic constant of lignin degrada¬ 
tion to char. K5: kinetic constant of lignin char degradation to volatiles. K6: kinetic 
constant of cellulose degradation to char. K7: kinetic constant of hemicellulose 
degradation to char. 



Fig. 9. Variation in percentage of the reaction kinetics between the samples 
between 26% and 30% of lignin. K2: kinetic constant of cellulose degradation to 
volatiles. K3: kinetic constant of lignin degradation to volatiles. K4: kinetic constant 
of lignin degradation to char. K5: kinetic constant of lignin char degradation to 
volatiles. K6: kinetic constant of cellulose degradation to char. 
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using the same kinetics parameters for the pure hemicellulose 547 

and cellulose and changing the parameters related with the lignin. 548 

This change implies that a relation between the three main compo- 549 

nents exit and justified the necessity of use a kinetic more complex 550 

than a first order kinetic to model the system. The used reaction 551 

pathway was the same as in pure samples. 552 

The sample treated with a mixture of ethanol/water (70-30) for 553 

1 h could not be adjusted by the same parameters. Consequently, a 554 

new set of kinetics parameters was needed to its fitting (average 555 

absolute deviation of 1.1%). Hence this experimentation is other 556 

proof of the structure of biomass has a role in its thermal degrada- 557 

tion. This deviation (A K = (/^No-Extracted — /^Extracted)//^Extracted) Can be 558 

observed in Fig. 7. In which it is presented that the extraction 559 

increases the degradation with temperature of the sample. The 560 

kinetic parameters are collected in Table 7. 561 

It is shown in Fig. 7 that the reactions more affected by the 562 

extraction are the degradation of the hemicellulose (Kl and I<7) 563 

and lignin (K3-K5). Cellulose (K2 and K6) decomposition is also 564 

increased but around 50% less. These results would be coherent 565 

with the expected behaviour due to the fact that lignin encloses 566 

cellulose and hemicellulose and it would be the most exposed. 567 

Regarding hemicellulose, it was expected that this process enhance 568 

its thermal breaking because it is water soluble. I<8 and I<9 (kinetic 569 

constant of cellulose and hemicellulose char degradation to vola- 570 

tiles respectively) were not affected. This result could be explained 571 

because the main contribution to the char comes from the lignin. It 572 

is also interesting the fact that the higher the degradation temper- 573 

ature is, the lower effect has the pre-treatment. Which would be 574 
expected because temperature enhances exponentially thermal 575 
degradation. 576 

4.2A.2. Grape skin. The main characteristic of this sample is that it 577 
does not have lignin. For this reason, the model had to reproduce 578 

the experimental behaviour using only the parameters related 579 

with hemicellulose and cellulose (Table 7). The model represents 580 

well the process (average absolute deviation of 1.9%) but there is 581 

a dramatic change of slope at 275 °C that cannot simulate. In addi- 582 

tion, it needed a change in the kinetics parameters which could be 583 

explained again by the effect of the biomass structure. 584 

4.2A.3. Hydrolysed samples. As we mentioned before, the grape 585 

seeds suffered for 1 h three different hydrolysis process at temper- 586 

atures of 250 °C, 300 °C and 340 °C. The result was a substance par- 587 

tially degraded with high content of charcoal. These three 588 

degraded samples were used in a TGA whose results were fitted 589 

with only a set of parameters (Table 7) with an average absolute 590 

deviation of 1.1% at 250 °C, 0.86% at 300 °C and 0.68% at 340 °C. 591 

This previous degradation generates a material which higher resis- 592 


0,16 

0,14 



Fig. 10. Kinetics constant in each isothermal process. Kl: kinetic constant of hemicellulose degradation to volatiles. K2: kinetic constant of cellulose degradation to volatiles. 
K3: kinetic constant of lignin degradation to volatiles. K4: kinetic constant of lignin degradation to char. I<6: kinetic constant of cellulose degradation to char. K7: kinetic 
constant of hemicellulose degradation to char. K9: kinetic constant of hemicellulose char degradation to volatiles. 
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tance against thermal degradation due to the fact that some char¬ 
coal was formed. This statement is shown in Fig. 8, where the var¬ 
iation is defined by the following mathematical expression 

K = (A/Civj 0 _Hydrolysed — ^Hydrolysed)/^Hydrolysed* If exhibits that the kinet- 

ics decrease respect the non-treated samples. However, hemicellu- 
lose kinetics presents the opposite behaviour but it is not 
representative because the main part of hemicellulose disappear 
during the hydrolysis. 

The lignin kinetic variations observed in Fig. 8 could be justified 
again with its protection function like in part 4.2.1.1. On the other 
hand, cellulose shows in this case a higher modification in kinetics 
due to their degradation during the hydrolysis. Finally, tempera¬ 
ture reduce again the differences between the thermal degradation 
of both samples. K8 and I<9 do not appear because of their low con¬ 
tribution to the total char. 

4.2.2. Woody biomass 

In this point, the thermal decomposition in isothermal and no 
isothermal conditions of woody biomass is considered. All the fit¬ 
tings needed a modification in the kinetics parameters. 

4.2.2.1. Non-isothermal process. Using the parameters obtained 
from the grape seeds the wood degradation could not be directly 
simulated because it has a different structure. The experimental 
data were divided in several groups depending on its lignin content 
(Table 4) to can use the same parameters for each group (Table 7). 
This division shows that a difference in the composition implies a 
change in the decomposition too. This could be explained by an 
interaction between the species in the wood, so when a greater 
amount of lignin exits the breaking of cellulose into volatiles and 
its charcoal formation are reduced and the degradation into vola¬ 
tiles of the later too. In addition, a higher production of volatiles 
from lignin and degradation into charcoal is achieved (Fig. 9, 
AI< = (I< w% - K 2 o%)/I< 2 o%y However, with a content of lignin 
greater than 30% the kinetics evolution changes. The kinetics of 
the cellulosic part of sample continues decreasing with the amount 
of lignin but the kinetics related to lignin decreases too. This 
behaviour could be caused by an increment in the thermal degra¬ 
dation resistance of the sample due to the fact that lignin, which is 
the strongest component against thermal degradation, would be 
more than 30%. 

It is remarkable that kinetic constants related with hemicellu¬ 
lose degradation do no change their values (for this reason they 
are not present in Fig. 9). This result could be explained by the fact 
that in this study case there were not previous factors that can 
solve it. So, its interactions with cellulose and lignin would be inde¬ 
pendent to the lignin concentration. The role of the temperature is 
the same as in parts 4.2.1.1 and 4.2.1.3 and K8 is not represented 
because there is no change in its value. 

4.2.22. Isothermal process. Finally, the TGA of a sample of Acer Sac- 
charum was studied in isothermal conditions at 150, 200, 250, 300 
and 350 °C. The adjustment needed a set of parameters for each 
temperature, when it is higher than 200 °C, and different from 
the parameters used in the non-isothermal process (Table 7). Their 
average absolute deviations were: 0.71%, 0.70%, 0.39%, 0.63% and 
2.68% respectively. This discrepancy between the kinetics due to 
the type of process could be caused by a protective interaction 
between species. This means that, in an isothermal mode, the 
decomposition is low up to a certain temperature is reached 
(250 °C) and hemicellulose degradation starts enhancing the 
decomposition of cellulose and lignin. Besides, there is an 
enhanced in lignin degradation 300 °C, when cellulose would start 
its degradation. This idea could explain the drastic change in the 
kinetics constants shown in Fig. 10 at 250 °C and 300 °C. 


It can be seen in Fig. 10 that the thermal degradation at isother- 655 
mal conditions depends on temperature strongly (as it was 656 
expected). In this case K5 and K8 are not present in the graph, 657 
because they did not change. This would be caused by the fact that 658 
the maximum operational temperature (350 °C) is not high enough 659 
to break lignin or cellulose char. 660 



T,ec 

Fig. Al. Fitting for the hemicellulose sample decomposition with a heating rate of 
20 °C/min. 



Fig. A2. Fitting for the cellulose sample decomposition with a heating rate of 
20 °C/min. 



Fig. A3. Fitting for the cellulose sample decomposition with a heating rate 
of 10 °C/min. 



Please cite this article in press as: Cabeza A et al. Autocatalytic kinetic model for thermogravimetric analysis and composition estimation of biomass and 
polymeric fractions. Fuel (2015), http://dx.doi.Org/10.1016/j.fuel.2015.01.048 


Fig. A4. Fitting for the cellulose sample decomposition with a heating rate 
of 5 °C/min. 
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Fig. A5. Fitting for the lignin sample from Turku decomposition with a heating rate 
of 20 °C/min. 



Fig. A6. Fitting for the alkaline lignin sample decomposition with a heating rate of 
20 °C/min. 



- XHC -TOTAL . EXP - XC - XL 

Fig. A7. Fitting for the grape seeds decomposition with a heating rate of 20 °C/min. 
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Fig. A8. Fitting for the extracted grape seeds decomposition with a heating rate of 
20 °C/min. 
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Fig. A9. Fitting for the grape seeds decomposition hydrolysed at 250 °C for 1 h with 
a heating rate of 20 °C/min. 
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Fig. A10. Fitting for the grape seeds decomposition hydrolysed at 300 °C for 1 h 
with a heating rate of 20 °C/min. 
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Fig. All. Fitting for the grape seeds decomposition hydrolysed at 300 °C for 1 h 
with a heating rate of 20 °C/min. 



Fig. A12. Fitting for the grape skins decomposition with a heating rate of 20 °C/min. 
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Fig. A13. Fitting for the lime decomposition with a heating rate of 20 °C/min. 
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Fig. A14. Fitting for the plane tree decomposition with a heating rate of 20 °C/min. 
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Fig. A15. Fitting for the holm oak decomposition with a heating rate of 20 °C/min. 
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Fig. A16. Fitting for the Arce Sadarium decomposition with a heating rate of 20 °C/ 
min. 
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Fig. A17. Fitting for the cedar decomposition with a heating rate of 20 °C/min. 
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Fig. A18. Fitting for the catalpa decomposition with a heating rate of 20 °C/min. 
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Fig. A19. Fitting for the almond tree decomposition with a heating rate of 20 °C/ 
min. 
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Fig. A20. Fitting for the Arce Sadarium isothermal decomposition at 150 °C. 
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Fig. A21. Fitting for the Arce Sadarium isothermal decomposition at 200 °C. 
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Fig. A22. Fitting for the Arce Sadarium isothermal decomposition at 250 °C. 


4.3. Composition estimation 

Once all the experimental data have been adjusted, the capabil¬ 
ity of the model to estimate the initial composition of the biomass 
was tested. The sample used to try this estimation was the TGA of 
non-extracted grape seeds. 

The prediction implies an optimization problem in which the 
difference between the experimental and simulated TGA must be 
minimized changing the values of the initial composition (12). 
The problem was limited by the following restraints. It was 
assumed that there is not initial charcoal in the sample and that 
the initial composition of water, oil, hemicellulose, cellulose and 
lignin were in the ranges showed in Table 5. The amount of inert 
compounds was obtained by balance to the total. 

min^|M Exp -M|j : ^jmin < m i < m jmax5J 6 (10) 

The calculated composition is not very accurate because in 
some compounds the deviation is high, for example the maximum 
deviation for water was 60.5% (Table 6). But taking into account 
that the values for the maximum and the minimum of each com¬ 
ponent were stablished in a general way, the prediction is good 
enough. In order to improve these values, more experimental com¬ 
positions would be needed to fix a better optimization range. It is 
interesting that the calculated cellulose composition is closer to 
the experimental one than the hemicellulose composition. This 
result could be caused by the fact that the oil vapourization and 
hemicellulose degradation can appear both between 250 °C and 
300 °C. 


5. Conclusions 
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Fig. A23. Fitting for the Arce Sadarium isothermal decomposition at 300 °C. 



Fig. A24. Fitting for the Arce Sadarium isothermal decomposition at 350 °C. 


An auto-catalytic kinetic model has been developed for thermo- 
gravimetric analysis with an average absolute deviation between 
the simulation and the experimental data lower than 7% in all 
the studied cases. This model can simulate the behaviour of very 
different samples (seeds, grape skins and trees) as isothermal pro¬ 
cess as non-isothermal process providing the composition profiles 
of their individual components too. In addition, the model can 
reproduce the effect of the heating rate in the decomposition using 
a non-Arrhenius’ dependence with the temperature. Due to the fact 
that the kinetics parameters change with the type of biomass it is 
deduced that the structure of biomass has a very important role in 
thermal degradation. Also important is the composition because 
some species can work as a shield that avoids the degradation of 
the others until their cleavage start. Finally, a preliminary compo¬ 
sition estimation were done, starting from a TGA curve and esti¬ 
mating the composition of the biomass material. This prediction 
has an acceptable accuracy especially for cellulose and lignin (dif¬ 
ferences lower than 7%). However, the prediction of the essential 
oil is trick and in order to increase model fidelity, more 
experiments would be needed, which would allow to stablish bet¬ 
ter limits for the optimization ranges of the initial composition and 
to improve the kinetics parameters. 
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Appendix A.l. Graphical adjustments 

A I. Pure samples 

AAA. Hemicellulose 
Fig. Al. 

A 1.2. Cellulose 
Figs. A2-A4. 

A 1.3. Lignin 
Figs. A5-A6. 

A.2. Complex samples 

A2A. Grape samples 
Figs. A7-A12. 

A.2.2. Woody biomass 

2.2.1. Non-isothermal process 

2.2.1.1. Content of lignin: 20% 

Fig. A13. 

2.2.1.2. Content of lignin: 26% 

Figs. A14 and A15. 

2.2.1.3. Content of lignin: 30% 

Fig. A16. 

2.2.1.4. Content of lignin: >30% 

Figs. A17-A19. 

2.2.2. Isothermal process 
Figs. A20-A24. 
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